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We propose and study the charge transport through single and double quantum point contacts 
setup between helical Majorana modes and an interacting helical Luttinger liquid. We show that 
the differential conductance decreases for stronger repulsive interactions and that the point contacts 
become insulating above a critical interaction strength. For a single point contact, the differential 
conductance as a function of bias voltage shows a series of peaks due to Andreev reflection of 
electrons in the Majorana modes. In the case of two point contacts, interference phenomena make 
the structure of the individual resonance peaks less universal and show modulations with different 
separation distance between the contacts. For small separation distance the overall features remain 
similar to the case of a single point contact. 

PACS numbers: 71.10.Pm,74.45.+c,05.30.Pr 


I. INTRODUCTION 

The recent discovery of topological insulators^ has 
spurred tremendous interest in the topological phases of 
condensed-matter systems. Topological systems in two- 
dimensional (2D) systems are characterized by their pe¬ 
culiar symmetry-protected gapless one-dimensional (ID) 
edge states in the presence of a gapped bulk^^ In time- 
reversal invariant (TRI) systems, two types of ID edge 
states are especially remarkable: 

On the one hand, helical Dirac fermions, whose spin 
is locked to the momentum, were first theoretically 
predicted^ and experimentally realized^ as the edge 
states of 2D topological insulators. Being rather insensi¬ 
tive to disorder, these edge states have promising appli¬ 
cations in the fields of nanoelectronics and spintronics. 

On the other hand, helical Majorana modes have been 
predicted to exist as the edge state of TRI topologi¬ 
cal superconductors^ ^ 18 ! 22 The current interest in the 
search for various Majorana modes^dH in condensed- 
matter systems mainly stems from their possible applica¬ 
tions in fault-tolerant quantum computing. While some 
experimental signatures for Majorana zero modes exist¬ 
ing as the end states of effective ID topological supercon¬ 
ductors have already been found^dJ, conclusive evidence 
in particular of their non-Abelian exchange properties is 
still actively sought for. 

In this paper we focus on the charge transport between 
a system of ID helical Dirac fermions and a system of he¬ 
lical Majorana modes, which are tunnel-coupled by one 
or several quantum point contacts. Due to the Coulomb 
interaction between the electrons, the low-energy proper¬ 
ties of the helical Dirac fermions is described by the he¬ 
lical Luttinger liquid theor y 19 ! 20 and is possibly realized 
in the InAs/GaSb experimental setup by Du’s grouped. 
On the other hand, the helical Majorana modes can 
to a good approximation be treated as free Majorana 


fermions. While strong interactions between the con¬ 
stituent electrons and holes^i may destabilize the Majo¬ 
rana modest the nearby superconductor screens mod¬ 
erate interactions effectively^ As long as the Majorana 
modes exist, they behave largely as chargeless particles 
and can be regarded as free. 

Similar tunneling phenomena in heterostructures 
have been discussed using renormalization group (RG) 
analysis^ - — and scattering formalis m 27 ! 30 ! 33 in the case 
of noninteracting lead(s). In this paper, we calculate the 
tunneling current by using perturbation theory in the 
coupling between the Majorana modes and the interact¬ 
ing helical lead, using the interacting helical lead Green’s 
function obtained by bosonization and the noninteracting 
Majorana Green’s function as the unperturbed propaga¬ 
tors. We use a scaling analysis to establish that the tun¬ 
neling term is the most renormalization-group relevant 
local perturbation in our system. We consider a finite- 
size topological superconductor with discrete helical Ma¬ 
jorana energy levels and assume for simplicity that the 
level separation is larger than the tunneling rate. With 
this assumption, we derive analytic results for the tunnel¬ 
ing current through one or two quantum point contacts, 
and obtain the current-voltage relation by evaluating the 
analytic results numerically. 

For a single quantum point contact with a noninter¬ 
acting lead, the tunneling current is the same (up to an 
extra factor of two in the differential conductance due 
to two spins) as for chiral Majorana fermions^ It shows 
periodic peak structures originating from the perfect An¬ 
dreev reflection in different Majorana energy levels. With 
increasing repulsive interaction strengths, corresponding 
to a smaller Luttinger parameter K , the differential con¬ 
ductance begins to decrease and eventually vanishes com¬ 
pletely at the resonance positions. This effect may partly 
explain why perfect Andreev reflection is difficult to ob¬ 
serve even despite the existence of Majorana zero modes 
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in quantum wire experiments — The quantum critical be¬ 
havior for the tunneling through a single point contact is 
similar to the charge transport with two helical Luttinger 
leads connected by a quantum do t 35 ’ 36 or the informa¬ 
tion leakage in the helical lead connected to a Majorana 
mode.— 

As in a real experiment, the tunneling may not be 
perfectly local, we shall also consider the effect of ex¬ 
tended point contact s—^ using a model involving two 
point contacts. For two nearby quantum point contacts, 
the distance between the two contacts determines the in¬ 
terference structures in the tunneling current . 53 i 54 For 
distances much smaller than the boundary length of the 
topological superconductor, the interference changes the 
shape of the individual peaks, giving features similar to 
Fano resonances with the overall magnitude^ remaining 
periodic. When the separation distance is comparable 
with the boundary length we see the overall magnitude 
also experiences some modulations related to the sepa¬ 
ration scale. Those interference features make the trans¬ 
port signature less universal and possibly modify the scal¬ 
ing behavior,— constituting another reason why perfect 
transmission is hard to observe. 

This paper is organized in the following way: in Sec. II 
we present the setup, the corresponding model Hamilto¬ 
nian, and the perturbation scheme. In Sec. Ill we use 
a scaling analysis to identify the tunneling term as the 
most relevant term in our system and make a comparison 
with other systems or different boundary conditions. In 
Sec. IV we present the main analytic and numerical re¬ 
sults for single and double quantum point contacts, and 
discuss their physical interpretations. We summarize our 
results and compare our approach to using scattering 
eigenstates in Sec. V. 


II. TUNNEL JUNCTIONS 
A. Proposed setup 

To realize both ID helical Luttinger liquids as well 
as ID helical Majorana modes, we adapt a proposal of 
Ref. [22j based on thin films of 3D topological insula¬ 
tors (3DTI) such as Bi 2 Se 3 or Bi 2 Te 3 , see Fig. [T] Helical 
Dirac fermions emerge at the sample edges due to the 
mixing of the top and bottom surface bands of the thin 
film . 42 ’ 43 Moreover, helical Majorana modes can be real¬ 
ized by sandwiching the thin film between conventional 
s-wave superconductors. For opposite signs of the super¬ 
conductor pairing functions on the top and bottom layer, 
and sufficiently strong proximity-induced pairing ampli¬ 
tude (greater than the mixing gap of the two surfaces), 
helical Majorana modes are indeed formed as the edge 
states of the thin film* 22 

In contrast to the converter between helical Dirac 
fermions and Majorana modes proposed in Refs. ism, 
we study the charge transport through tunneling junc¬ 
tion^) between a helical Luttinger liquid and helical Ma- 



FIG. 1: Proposed setup for realizing heterojunction(s) of heli¬ 
cal Luttinger lead and helical Majorana modes using thin film 
of three dimensional topological insulator. The top and bot¬ 
tom superconductors’ order parameters carry different signs 
as proposed in the Ref. I22I . 



FIG. 2: Schematic figure of single tunneling junction 


jorana modes connected to ground. For systems with bro¬ 
ken time-reversal symmetry, where the helical edge states 
are replaced by chiral ones, such transport phenomena 
were studied with noninteracting lcads^H and interact¬ 
ing leads.— The schematic diagram for single tunneling 
junction is depicted in Fig. [2j where the voltage differ¬ 
ence between the two leads is controlled by the chemical 
potential /1 imposed on the helical Luttinger liquid lead. 
The tunneling amplitude t is controlled by the width of 
the junction and is related to the wavefunction overlap 
between the two leads. 

In the real experiments it may not be easy to fix the 
relative phases of two adjacent superconductors nor fine 
tune the chemical potential to the topological regimes. 
For the heterostructure setup we can make the helical 
Majorana modes by different types of realizations^, or 
change the helical modes to time reversal preserved dou¬ 
ble Majorana end states^. The results for tunnel trans¬ 
port to double Majorana end states is the limiting case of 
helical Majora mode with energy level separations going 
to infinity, as is shown in Fig. [5] 
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B. Model Hamiltonian 


We consider one helical Luttinger liquid lead and one 
helical Majorana fermion lead. The Hamiltonian describ¬ 
ing this system is H = H L + H Mo + )T) Q HT a + $H. The 
Hamiltonian Hr, — f°° dx'Hr, for the helical fermions is 
a Luttinger liquid Hamiltonian density, 


Hl — ivF 


ipl(x)d x ijj L (x) - il’Ux)d x ip R (x)} (1) 


-H(x) + ^ R {x)ip R (x) 

+U2ip' L {x)i>L {x)i’h {x)i>R ( x ) 


TJ" WVV (x)lpt (x)lp r (x) 

r=R,L 


The Hamiltonian for the grounded propagating Majorana 
fermions on a ring of circumference L is 

Hm 0 = j Y, / dx (vM,<j'y<r(x)d x j< T (x)) (2) 

cj J o 

Here vm ,,t = sgn(er)uM- The single particle tunneling 
term between the helical Luttinger liquid lead and helical 
Majorana fermion lead is described by^ 

H T=iYl [&-CTQ:Vv(Za) + (3) 

r,«r,a v 

Here, t rc m is the tunneling strength, and £ raa are com¬ 
plex numbers with Ifro-al = 1. r indicates the left/right 
movers in Hr, a denotes the spin index of the Majorana 
fermions, and a = 1 ,..,N is the number of tunneling 
channels (junctions) and x a /y a are their spatial coor¬ 
dinates in Luttinger/Majorana leads. We restrict our 
discussions to N = 1 and N = 2 in this paper but the 
extension to arbitrary N is straightforward and similar 
to the N = 2 case. The N = 2 case is illustrated in 
Fig. |3l 

The remaining SH term contains the leading 
instabilitie s 27 ! 28 under the renormalization group anal¬ 
ysis in the low-energy sector. We show in the next sec¬ 
tion why they are not important in our setup. With 
this simplification the full low-energy effective Hamilto¬ 
nian becomes H ~ H L + H Mo + H Toi , describing the 
single-particle tunneling between spinful Luttinger liq¬ 
uids and Majorana fermions lead. The tunneling charge 
current through site x a (in the helical Luttinger liquid 
coordinate system) is obtained by 



FIG. 3: Schematic figure for two tunneling junctions sepa¬ 
rated by spatial distance Ax and Ay 


are less than the coherence length). We choose a time 
dependent gauge transformation to move the chemical 
potentials in Hl to Hr by writing i^r/l —> e l,lt ipR/ L . 
By defining the Keldysh contour ordered Green’s func¬ 
tion G„ tR / Lt0l {t,t') = -i{T c {') <J (y a ,t)il)' R/L (x a ,t')}) we 
rewrite the particle current as 

m/e = St[ J2 t jaa e-^G< Jta (t,t)]. (4) 

j=R,L;<j 


This lesser mixed Green’s function Gf j (t, t) is ob¬ 
tained by perturbation theory as 


/_»y+ 1 r r 

G„, R/L , a (t,t) = J2 dT i--- M T c{ 

i=o L Jc Jc 

Tt (.Vot > tyHint {t\ ) ■ • ■ Hint (l~l')'4 ! R/L (Xa ; t )}) (5) 


In applying the Wick theorem in the Eq. (J5J) we should 
also include all possible four fermions interactions term 
(u 2 and U 4 term in the edge states Hamiltonian) be¬ 
tween any two fermions operators. We use the spinless 
bosonization^i as a way to sum up all orders of pertur¬ 
bations in the four fermions interactions on the Keldysh 
contour. The edge state correlators evaluated this way 
is thus fully dressed by the four fermions interactions in 
our treatment and we do not specify this aspect in the 
expression of Eq. ©. 

We bosonize the helical Luttinger liquids lead opera¬ 
tors by writing the fermion fields as: 


<4J 


* e (E V40)</v( x),H]} 

r 

~ e (Y (€r* a i>r(x a ) ~ Caa^liXa))) 

r,<j * 


^r{x) 

$l{x) 


^/2nao 

1 

y/Zirao 


V Re 


-i\Z4K4>R(x) 


y l& 




( 6 ) 


The total tunneling current is the coherent sum of 
the current from all tunneling channels (under the as¬ 
sumption that the separation distance between junctions 


with r\RiR as the Klein factor chosen to satisfy the 
fermion anti-commutation rule and ao as the lattice spac¬ 
ing cutoff for the linear spectrum. We define the bosonic 
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fields <J>, 0 = </>l ± (j)R and rewrite H 0 = Hl + Hm 0 and 
H t as 

v r°° i 

Ho = ^J dx:[K(d x Q) 2 + -(d x $) 2 }: 

+ / dxj a (x)d x "f a (x) 

a J 0 

Ht = 

era. v 

+ (7) 

with Luttinger parameter K = , anc j veloc- 

° ^ Y ZlXVF-\-U4.-\-U2 

ity v = v F y/(l + 2^:) 2 - (^j) 2 . Eq. Q serves as the 
main Hamiltonian for computing the tunneling current in 
section llVl For a single tunneling point contact with time 
reversal symmetry preserved we set t R -\ a = fz 4 a = t a 
and otherwise zero. We discuss why other relevant per¬ 
turbations 5H are not important in this time reversal 
preserved system in the next section. 


III. SCALING ANALYSIS 


Following the discussions in Ref. [2?] for a single tun¬ 
neling junction located at x = 0 , the most relevant terms 
SH other than the tunneling term Ht are 

SH = Vi[4(0)^ fl (0) + 4(0)^i(0)] ( 8 ) 

+ [E 2 ^( 0 )t/>L( 0 ) + Ai/) R (0)i/jl(0) + h.c.]. 


Here, the Vj terms represent the chemical potential 
change due to the presence of the tunneling junction 
(also called quantum point contact). The V 2 terms stand 
for backscattering due to the point contact and A is the 
Cooper pair gap magnitude induced at x = 0 via prox¬ 
imity effectj^iSS Rewriting the fermionic operators via 
Eq. © we get 6H in bosonized form as 


SH = -^$(0) 



+ 



(9) 


with A = |A|e*A The V\ terms can be absorbed in the 
definition of $( 2 :) by the shift $(x) —¥ <h(x) — 2 ^^ sgn(x). 
For the rest of the terms in Ht and SH the scaling di¬ 
mensions around Hq are^ 



D[V 2 \ = K, 

D[ A] = l/K. (10) 

The term i in D[t\ comes from the scaling dimension 
of helical Majorana modes D[ 7 ] = i, which is the same 


as chiral ones, assuming its spectrum is continuous (or 
the boundary of the topological superconductor being 
infinite). For a time reversal symmetric Hamiltonian 
the backscattering term proportional to V 2 , being the 
only relevant term in the repulsive interaction regime 
(0 < K < 1), is forbidden. For repulsive interactions 
1 < D[t\ < D[A] and thus the most important terms 
(marginally relevant) in perturbation for Ht+SH around 
Hq is the tunneling term Ht- 

For the short topological superconductors considered 
in this paper, the helical Majorana edge states become 
discretized and £>[ 7 ] ~ 0. Under this approximation 
D[t\ ~ \ (K + -^) becomes relevant for 2— \/3 < K < 1, 
while -D[A] stays irrelevant in the repulsive regime, in¬ 
dicating the same quantum phase transition (metallic to 
insulating) as for helical Luttinger liquids connected via 
a quantum dot ^ 6 in the repulsive regime. 

For different geometries, such as a Luttinger liquid ter¬ 
minated at a Majorana zero mode end stated or helical 
Luttinger liquid connected to a time-reversal breaking 
topological superconductor (with chiral Majorana modes 
as its edge state) ^ the backscattering V 2 term is rele¬ 
vant for K < 1 and the low-energy physics is determined 
by a new fixed point Hamiltonians^ 

h' 0 =H 0 -sin(V4^E>(0)) , ( 11 ) 

7rao \ / 


which fixes the value of <h(0) = for V 2 > 0. Under 

this constraint the scaling dimension of tunneling term 
D [t] with D [ 7 ] ~ 0 becomes D [t] ~ ^ and is relevant for 
1/2 < K < 1, giving rise to the transition between per¬ 
fect normal and perfect Andreev reflection at K = 1/2 in 
this system. 27 - The transition from perfect normal to per¬ 
fect Andreev reflection is shown as insulating to metallic 
transition in the charge transport. The key difference 
from our setup is the different scaling behavior (different 
power law dependence), controlled by the density-density 
interaction strength in the helical Luttinger liquid, in the 
differential conductance as a function of bias voltage or 
temperature. 


IV. EVALUATING THE CURRENT 

In this section we carry out the calculation of the cur¬ 
rent for a single point contact and a double point con¬ 
tact at zero temperature. We start by finding analytic 
expressions for the helical Luttinger liquid and dressed 
helical Majorana modes Green’s functions. From there 
we compute the current numerically by using Eqs. Q and 
©, and thus obtain the current-voltage relation numer¬ 
ically. By taking the derivative numerically we get the 
differential conductance as a function of voltage. We find 
a metallic to insulating quantum phase transition (near 
zero bias) with increasing repulsive interaction, and less 
universal patterns owing to the interference nature in the 
case of double point contacts. 










5 


As a side remark, notice that the computation car¬ 
ried out here is not the one loop RG calculations men¬ 
tioned in the previous section. We perform a diagramatic 
based resummation of perturbative terms and the eval¬ 
uated differential conductivities depend explicitly on the 
choice of linear momentum cutoff A. The choice of A 
depends on the particular realizations of the helical Lut- 
tinger modes, i.e. material dependent. The cutoff de¬ 
pendence, as shown in Appendix O is consistent with 
the trend we expect from usual higher order (two loops 
or more) of RG calculations. That is, for larger cutoff 
A, the deviations from what we expect from lowest RG 
analysis are larger. In the rest of the paper we choose 
A = 10~ 2 e_F as a typical value of modeling the lineariza¬ 
tion of some quadratic bands at the Fermi surface, or 
the band touching point where the edge states of 2DTI 
become mixed with the bulk band. 


A. Single point contact 


We start with single point contact between the heli¬ 
cal Luttinger liquid and helical Majorana modes realized 
in a time-reversal symmetric topological superconductor. 
From Ref. [H the Keldysh component of bare (uncou¬ 
pled) lead Green functions, defined as G^ l/r {t,t') = 

—i{T c {i/ L / R {r)il) lf L ^ R {T')e~' LllL / R ^ T '~ T ' ) }), expressed in fre¬ 
quency space at zero temperature are^l 


n++ 

L/R 


M = 


x 

g V l/ M = 

X 


G E» = 




-mI 


2k-1 


47t 2 u 2k F(2k) 

(Ji(k)9(uj — (i) — h(n)9(g — ui)^J 


r ( K ) 2 


-Ml 2 


( 12 ) 


47t 2 v 2k T(2k) 

h*{n)9{u — /j.) + h*(n)9(g, — w)^ 

v 2K I (2k) 


v 2k T(2k) 1 


- - n) 


Here k = ^(AT+l/AT) and Ji(k) = 2e mK sin(7TK)r(l — k) 2 
and plus/minus sign on G^, L/R indicates its labeling on 
the Keldysh contour (with G^ as time ordered and 
G 4>l/r as an tiAime ordered). By relabeling the spin index 
in 7 ^ by the left/right-movers label of the Luttinger lead 
operator, the steady state charge current is expressed as 


(I) = e3? 


E 


t t* 
h J, nU j,m 


2 


+ G^ j nm {u})G^. (w)) 



Here n, m denote the discrete energy levels in the finite 
size helical Majorana modes. Ea. (fT3l) follows from main¬ 
taining the structure of first order expansion in Eq. ([5|) 


and resum all higher order terms through the ’’dressed” 
helical Majorana Green’s function. The retarded helical 
Majorana Green’s function contains higher order terms 
through inclusion of self energy terms: 

G« in Ju>) = G^» (14) 

+ E g E>)E«'M g *,,>)- 

Here the ’’bare” retarded helical Majorana Green’s 
function is G^^(w) = 5 n , m /{w — e n ,j + *0 + ) with 
e n ,j = hv M sgn(j)^ (with sgn(j) = +/- for 

L/R), and the retarded self energy is T,^ nm {ui) = 

qr ^ given by the Dyson equation. The 
dressed helical Majorana lesser Green’s function 

is G 7 3 -,n m M = G Ei( w )E«'H G tcJ a; ) with 

G< ( uj ). Similar expressions hold for 

G Em( w ) and G r 

In Eq. (1131) the summation over integers n, m (and 
1,1' in Ea. dTll) ') refers to the sum over discrete Majorana 
modes energy level indices. From this the expression 
for the current is related to the evaluation of G^.(u;) = 

G J“M, G/p. (w) = G^ + (w) - G" + (w), and the afore¬ 
mentioned dressed helical Majorana Green’s functions. 
In this paper we assume the energy difference between 
different Majorana modes is sufficient large (greater than 
the broadening effect coming from coupling with the Lut¬ 
tinger lead) such that Eynm(w) ~ T,j{Lo)5 n ^m to simplify 
the calculation. In other words, we consider the helical 
topological insulator as short, such that the finite size 
makes the energy difference between discrete Majorana 
modes sufficient large so that the overlap between them 
is negligible. Under this assumption the helical Majo¬ 
rana Green’s functions G> nm (w) and G ^ (w) are diag¬ 
onal and the analytic expression for Eq. m is obtained. 
From there we evaluate the current numerically and ob¬ 
tain its relation with bias voltage V = (/x — 0) /e to eval¬ 
uate the differential conductance. The results for single 
tunneling junction, with different Luttinger parameters 
K indicating different interaction strengths, are shown in 

Fig- SI 

In Fig. SI we see the perfect transmission (maximum 
differential conductance) at zero voltage for a noninter¬ 
acting (K = 1) helical Luttinger liquid. It originates 
from the perfect Andreev reflection between the metal 
and superconductor mediated by the helical Majorana 
modes.— The Majorana modes inside the superconduct¬ 
ing gap serve as resonance levels which facilitates the 
Andreev reflection process and give a differential conduc¬ 
tance value cqe 2 /h with g\ = 2x2 = 4, reflecting particle- 
hole and spin symmetry. This perfect transmission sig¬ 
nature is used to identify the Majorana zero modes in 
the nanowire experiments.— The periodic peaks at finite 
bias voltages, similar to the case of tunneling measure¬ 
ment chiral Majorana modes discussed in Ref. [30], come 
from discrete Majorana energy levels with energy differ- 
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dI/dV(gi^ 2 /h) 



eV/A 


dI/dV(g,e 2 /h) 



FIG. 4: Differential conductance as a function of voltage for 
different Luttinger parameters in the helical Luttinger lead 
connected by single quantum point contact with grounded he¬ 
lical Majorana modes. The Luttinger parameters are: K = 1 
(blue solid), K = 0.5 (purple dot dashed), K = 0.3 (brown 
dashed). Other parameters are t = 0.05A, the length of the 
edge of topological superconductor L = 10 3 ao, and linear 
spectrum cutoff A = 10 _2 £f = 10 

ence (peak intervals) set by the physical size of the edge 
of helical topological superconductor. For (repulsive) in¬ 
teracting lead the general feature is the suppression of 
the resonant conductance peaks, both in peak magni¬ 
tude and width, and the spreading out of spectral weight 
away from the resonance levels/ 4 '' 

The spreading of spectral weight makes the transition 
from perfect Andreev reflection to perfect normal reflec¬ 
tion more difficult to observe at smaller Luttinger pa¬ 
rameter K. This is because the Majorana levels begins 
to merge together (as shown for K = 0.3 case in the 
Fig. SI), violating our starting assumption that the levels 
are sufficient far apart. To illustrate this kind of metallic 
(perfect Andreev reflection) to insulating (perfect nor¬ 
mal reflection) behavior in this single tunneling junction 
context, we plot the differential conductance for a single 
level (zero energy) helical Majorana mode in the Fig. [5] 
We see that the transition takes place between K = 0.3 
and K = 0.2 with marked tendency differences between 
the two at finite bias. 

The scaling analysis mentioned in the previous section 
for single Majorana mode gives D[t\ = {K + l/AT)/4, 
resulting in a critical Luttinger parameter K cr = 2 — 
\/3 for the repulsive helical Luttinger lead. This scal¬ 
ing/criticality behavior is the same as the case for two he¬ 
lical Luttinger liquids connected (with particle-hole sym¬ 
metry imposed, or = -~fi 2 = eV/2) by a noninteract¬ 
ing single level quantum dot discussed in Ref. [36j]. The 
vanishing charge transport at zero bias below K = K cr 
corresponds to the transition point where some quantum 
information stored by qubits formed by Majorana modes 
is maintained and does not decohere completely— In gen¬ 
eral, for charge transport we can formally make analogy 
between the two helical Luttinger leads with particle-hole 
symmetric driven voltage connected via a noninteract¬ 
ing multi-level quantum dot system with our single heli¬ 


FIG. 5: Differential conductance as a function of voltage for 
different Luttinger parameters with a single Majorana level: 
K = 1 (blue solid), K = 0.5 (purple dot dashed), K = 0.3 
(brown dashed),A = 0.2 (green dotted). Other parameters 
are the same as in Fig. [4] 

cal Luttinger lead connected with the helical Majorana 
modes. 

Note that once introducing time-reversal-invariance 
breaking terms^ or different ways^i of connecting the 
Majorana modes, the critical behavior in the charge 
transport could occur at different K cr . For example, if 
we replaced the helical Majorana modes by the chiral 
Majorana modes (the edge state of a topological super¬ 
conductor with broken time reversal symmetry)^ simi¬ 
lar metallic to insulating behavior is seen but with scal¬ 
ing behavior controlled by D[t\ = 1/2AT, or K cr = 1/2. 
Thus for single Majorana modes (or other modes suffi¬ 
ciently far apart such that the overlap is not significant) 
different scaling behaviors in the single tunnel junction 
transport reveals a great deal of information about the 
boundary conditions imposed on the Luttinger liquid. 

In a real experiment there is a finite length region 
where tunneling between the Majorana modes and helical 
lead occurs. Assuming a separable form for the spatial 
dependence of tunneling ter m 39 : 40 the analytic expres¬ 
sions for weak tunneling current are obtained for the case 
of leads made of same type of material. For tunneling be¬ 
tween helical Luttinger leads, the power-law dependence 
is modified for an extended contact compared to the case 
of a point like contact.— Extension of this formulation 
to infinite-size helical Majorana lead case seems to be 
straightforward but not so easy for the case of a finite 
helical Majorana lead. Therefore, we proceed with the 
simpler case: the case of two quantum point contacts. 

B. Double point contact 

Proposals for double point contact setups are mainly 
related to the study of quantum interference effect s 44 : 45 
and the quasi-particles statistics of the edge states— 
Similar setups for the helical edge states have been 
discussed^ 52 » 54 with the possible applications for elec¬ 
tronic means of spin pumpin g 53 : 55 . Different types 
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of interferometers have also been proposed for het¬ 
erostructures of topological superconductor and normal 
metal/edge states of topological insulators such as the 
Majorana Dirac converter! 22,23 ! 49 Here, the double point 
contact between the helical Luttinger liquid and the Ma¬ 
jorana modes is yet another type of heterostructure show¬ 
ing the quantum interference which is analogous to the 
two point source interference in optics. 

For two point contacts the total current passing 
through those contacts is (I) = ( I X1 + I X2 ). Without 
loss of generality we choose X\ = 0, X 2 = x, jq = 0 and 
y 2 = y■ Here Xi denotes spatial coordinate of the fermion 


operators in the helical Luttinger lead and yi denotes that 
of the helical Majorana operators in the topological su¬ 
perconductor. We evaluate the current (I) via perturba¬ 
tions on the tunneling term Ht on the Keldysh contour. 
To simplify the notation we denote ^ as the bare 

(unperturbed) helical fermion/Majorana mode with chi¬ 
ral (or spin) and position index a and concentrate on the 
structure of perturbation in the Dyson equation without 
bookkeeping the Keldysh contour labels on the Green’s 
functions for the moment. We get: 


G la = g£> + G^\t a \ * 2 G^G la + G™t a + G^tpGwXG^ + GW Jfi M 2 G (15) 

G lalp = G% p +G^t a G Ml ^G lf) + G^ 7 p \tp\ 2 G^G lft + G^ 7p tpG^J* a G^ + G^fG^G^W) 


Note that we do not have spin flip process in the tun¬ 
neling term, and the G^^ a or G 7/37q , are diagonal in spin 
space and functions of differences in the spatial coordi¬ 
nate. In these simplified notations the current at position 
x a of the helical Luttinger lead coordinate is 

(t„) = e3 ?[J du> (| t a \ 2 G la Gip a + t a tpG 7a7p G^,^ a ^\17) 

I 


It is easy to check that above formula gives the single 
point contact result (ThU) by taking t a = t/y/2, tp = 0, 
and with the Langreth theorem^ (to denote the contour 
order). Following the same recipes the current for two 
point contacts with t a = and tp = t 2 /\/2 is then 

expressed as: 


(-0 — + {-^ 2 , j )} 


j=± 1 


(G 7 .(uj)G^.(uj) + G 7 .(uj)G^. (w)) + t!t* 2 (G 7j (u, -yi 2 )G^.(u],xi 2 ) + G 7j (u, -y\ 2 )G^.(uj,xi 2 
(1x2, j) = o/ du(\i 2 \ 2 (G*(u;)G<.(u) + G<. (w)G$. (w)) + t 2 tl(G^. (w, 2/12 )G^.(u;, -x\ 2 ) + G<.(uj,yi 2 )G$,.(u], -x 12 


r 


Her eG^Vz/) = E *" Sgn0) 


and G^ )< (w,y) = 




2iri e~ lsgn ^' >2j ^ ±v d(—.Lo)S(e n j — w) are the unper¬ 
turbed retarded and lesser Green’s function for helical 
Majorana modes, label j = ±1 denotes left/right moving 
mode, and x \ 2 = X\ — x 2 and y\ 2 = y\ — y 2 are the spatial 
coordinate differences. Following Eqs. (fl5l) and (fTBl) and 
the Langreth rule we obtain the various dressed Majo¬ 
rana Green’s functions and unperturbed helical fermions 
Green’s functions needed for evaluating the current. The 
derivations and analytic expressions for various Green’s 
functions are shown in the Appendix [5] and [C] With the 
analytic expressions shown in the Appendices, we per¬ 
form numerical integrals to compute the current (1181) and 
obtain the differential conductance by taking numerical 
derivatives with respect to the source drain voltage V. 


The results are shown in the Figs. [HI [HI 

We chose a small separation length (xi 2 = y\ 2 = 
10 _2 L) between the two contacts in the Fig. H3 and fixed 
the tunneling strengths of the two point contacts to 
be identical. The interference effect due to two point 
contacts for the weakly interacting lead (K ~ 1) is 
not apparent, and the resonance structure is similar 
to the single point contact. For a noninteracting lead 
(K = 1 or blue solid line in Fig. [6J the differential con¬ 
ductance reaches its maximum value g 2 e 2 /h with g 2 = 

2 (particle hole) x 2 (spin) x 2(2 tunneling points) = 8 
when the chemical potential of the helical lead is in line 
with the discrete Majorana energy levels. For a heli¬ 
cal lead with stronger repulsion (say K = 0.5 or purple 
dashed line in Fig. [Gj) we see features similar to the single 
point contact (with shrinking peak width and height at 








dI/dVfe 2C 2 /h) 



dI/dVfe 2C 2 /h) 



FIG. 6 : Differential conductance v.s. voltage for different 
Luttinger parameters K = 1, K = 0.5, and K = 0.3 with 
separation distance 212 = 3/12 = 10ao. We choose the tunnel¬ 
ing term ti/A = fa/ A = 0.01, and the length of the helical 
Majorana modes L = 10 3 ao with ao denoting lattice spacing. 
A = 10 ' 2 £f = lOhvM is the linear spectrum cutoff in the 
helical Luttinger lead. 

dI/dV( S 2 c 2 /h) 



FIG. 7: Differential conductance v.s. voltage for Luttinger pa¬ 
rameter K = 0.5 with different separation lengths. We choose 
*12 = 3/12 with £12 = 3.18ao (blue solid), £12 = 7.96ao (purple 
dot dashed), and 212 = 10ao (brown dashed). Other param¬ 
eters: fi/A = fa/A = 0.01, L = 10 3 ao, and A = 10 _2 £f = 
10 hvM^r-- Top right inset shows the enlarged figure for peaks 
around eV/A = 0.3. 


resonance value and transfer of spectral weight away from 
the resonance) and the effect of interference between two 
point contacts. Around zero bias the peak splits into two, 
similar to the physics of Fano resonance, and slight mod¬ 
ulations in the resonance positions in other finite voltage 
peaks. To further study the interference effect we fix 
K — 0.5 and plot different separation lengths (still keep¬ 
ing X 12 = yi 2 and X 12 /L ~ 10 -2 ) in Fig. [7] We see the 
subpeak structure (see the inset of Fig. □ also emerges 
nearby finite voltage resonance peaks with peak heights 
at a fixed voltage depending on the separation distance. 
This kind of subpeak structure mainly comes from the 
change in the real part of self energy correction on the 
Majorana Green’s function, which emerges with the can- 


FIG. 8 : Differential conductance v.s. voltage for different 
Luttinger parameters K = 1 (blue solid), K = 0.5 (purple 
dot dashed), and K = 0.3 (brown dashed) with separation 
distance 212 = 1/12 = 10 2 ao- We choose the tunneling term 
ti/A = t 2 /A sb 0.01, and the length of the helical Majorana 
modes L = 10 3 ao with a 0 denoting lattice spacing. A = 
10 ~ 2 £f = 10 hvM is the linear spectrum cutoff in the helical 
Luttinger lead. 


dI/dV(g 2 i? 2 /h) 



FIG. 9: Differential conductance v.s. tunneling amplitude 
phase difference (j> for different Luttinger parameters K = 1 
(blue solid), K = 0.5 (purple dot dashed), and K = 0.3 
(brown dashed), ti = t 2 e 1 ^ with |ti|/A = 0.01. Other pa¬ 
rameters: 212 = yi 2 = 10ao, A = 10~ 2 cf = 10and 
L = 10 3 a 0 . 


cellation of fast oscillating term related to e lkpsgn ^ Xl2 in 
the G in the Appendix [Cl from different orientations. 

For longer separation distance the interference effect 
also brings a change in the peak heights. To demon¬ 
strate this we choose x ±2 = 3/12 = 0.1L in Fig. [ 8 ] such 
that the separation distance is one tenth of the linear 
dimension of the Majorana modes. Other than the sub¬ 
peak structure seen for K = 0.5 case we now also see 
modulations in the resonance peak heights. This larger 
envelope (modulation with large voltage range) is associ¬ 
ated with the separation length scale being comparable 
with the Majorana system size. For off resonance region 
(K = 1 plot with voltage between eV /A = 0 to 0.1 in 
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Fig. m for example) we also see small oscillations around 
zero which is attributed to the inaccuracy of numerical 
integrals for fast oscillating functions. 

For general double point contacts we could have differ¬ 
ent tunneling amplitudes 1 i, t 2 and different separation 
distances \x\ 2 \ ^ I 2 / 12 1- For |H| » \t 2 \ or |ti| < |t 2 | the 
transmitted current is dominated by one of the point con¬ 
tact, and the result is basically the same as that of the 
single point contact. For |xi 2 | ^ \yi 2 \ but with |fi| = |t 2 | 
and |a:i 2 1 ~ I 3/12 1 the general features are similar to what 
we have mentioned in this section. Here we discuss the 
case of identical separation length |a:i 2 1 = I 2 / 12 I = 10ao 
but with different tunneling amplitude t\ = t 2 e x<l> . We 
plot differential conductance around zero bias as a func¬ 
tion of the tunneling phase difference </> for different Lut- 
tinger parameters in Fig. [9] For small separation distance 
chosen here </> = 7rorti=— 1 2 leads to almost complete 
cancellation of the resonance peak. For larger separa¬ 
tion distance the general feature is the same (decreasing 
dl/dV with increasing <J)) but with finite conductance 
even at <f> = it. 


V. CONCLUSION 

We have investigated the charge transport between a 
helical Luttinger liquid and a system of helical Majo- 
rana fermions coupled by single and double quantum 
point contacts. The helical Luttinger liquid is realized 
as the one dimensional edge state of a thin film of a 
3D topological insulator with the inclusion of short-range 
repulsion. The helical Majorana fermion could be real¬ 
ized in noncentrosymmetric topological superconductor 
or proximity-induced effective topological superconduc¬ 
tor with time reversal symmetry. For a single tunneling 
point contact we find that perfect Andreev reflection oc¬ 
curs only for a noninteracting helical lead. Increasing the 
repulsive interaction strength leads to the suppression of 
the differential conductance on resonance and shifts the 
weight away from resonance. This feature is similar to 
the case of two Luttinger leads connected by a noninter¬ 
acting quantum dol— with particle hole-symmetric bias 
voltage (fj, 1 = —g 2 = eV/2). 

We then studied the case of two quantum point con¬ 
tacts. For small separation distance (x\ 2 -C L with L 
being the size of the edge of the topological supercon¬ 
ductor), the interference from the two point contacts 
strongly changes the shape of the individual resonance 
peaks but does not affect the overall magnitude at dif¬ 
ferent Majorana mode energies. At larger separation dis¬ 
tance (X 12 ~ 10 _1 L) we observe modulations in the mag¬ 
nitude and shape of individual resonance peaks resulting 
from two point interference. 

In a real experimental setup, the point contact may 
not be perfect, in which case an extended contact may 
provide a better description ; 39 ' 40 The analytic results of 
the perturbation theory in the tunneling get more com¬ 
plicated with an increased number of tunneling channels, 


as shown for the case of two point contacts in this paper. 
We conjecture, based on our result at small separation 
distance, that with sufficiently small size of this extended 
point contact (Ax ~ Ay -C L), the overall transport be¬ 
havior will be similar to the single point contact. The 
detailed scaling behavior— or the shape of the individual 
resonance peak can be different and the transport signa¬ 
ture gets modified by the interaction more significantly. 
This can also be viewed as a generalization of the scaling 
behavior change due to the modification of the boundary 
conditions as mentioned in Ref. [28], 

As a final remark, the noninteracting limit (K = 1) 
of our results can also be derived by the scattering func¬ 
tion formalism as done for the chiral Majorana case— 
For repulsive interactions (K < 1), one can in principle 
use the Bethe ansatz scattering eigenstate s ^ ^ and de¬ 
rive the tunneling current for a single point contact. This 
formulation might be an extension of the perturbative ap¬ 
proach introduced here, had the issues of complex Bethe 
momenta be clarifiedi 5 ' Different type of interacting leads 
realizing different kind of Luttinger liquids^ can also be 
connected with noninteracting Majorana modes, which 
leaves unique transport signature due to different scaling 
behavior for an ideal single point contact. 
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Appendix A: Cutoff dependence of differential 
conductance 

In this section we evaluate numerically the zero 
bias differential conductance at various different cutoffs, 
maintaining the cutoff energies at the order of 10 _2 £f- In 
the main text we choose A = 10 (i.e. with ep/hvF = 10 3 
inverse length unit) and tunneling amplitude t = 0.5 in 
Fig. S]and Fig. [5] and different t for the rest of the fig¬ 
ures. Here, for demonstration purpose, we fix t = 0.5 
and all other parameters the same as those in Fig. 0] 
and Fig. [5J and we vary A from 10 to 40 stepped by 
10. The result is shown in log-log scale in the Fig. [TO] 
For Kcritical — 0.26 < K < 1 the lowest order RG 
flow to the fixed point of perfect transmission, giving 
dl/dV |y_>.o = 1 (in unit of g\e 2 /h). 

In Fig. E] the calculated zero bias differential conduc¬ 
tance depends explicitly on the linear momentum cutoff 
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dI/dV | v=0 



FIG. 10: Zero bias differential conductance v.s. linear mo¬ 
mentum cutoff A for fixed tunneling amplitude t. 


A, with larger A giving smaller value and a generic trend 
of decreasing transmission amplitude for decreasing K 
(or stronger repulsion in helical Luttinger lead, shown for 
K = 0.7 to 0.3). For generic two loops RG (or higher or¬ 
der) the next leading correction normally takes the form 
of (1 + 1/F(A)), with 1/F(A) -> 0 as A -> 0. The ex¬ 
plicit form of the higher order corrections F( A) depends 
on the specific Hamiltonian. The trend we see in Fig. [TO] 
is consistent with the naive higher order RG. However we 
shall bear in mind that this is not the RG type of calcu¬ 
lations, but a fixed cutoff with inclusion of most of the 
perturbative terms (neglecting the level crossing terms) 
via the Dyson’s approach. 


Appendix B: Luttinger lead correlators 


For both t 2 and t\ on the top branch, or time ordered 
branch, we get 


F (1)++ (f,:r) 


F {2)++ (t,x) 


9(t)F^- + {t,x) + e(-t)F^ + -(t,x) 


-1 


In 


x 2 + (a + iv\t\) 2 


27r 

e(t)F^~+(t,x) 


(B5) 


9(-t)F^+-(t,x) 


-1 

2ir 


■In 


a + iv\t\ — isgn[t]x 
a + iv\t\ + isgn[t\x 


(B 6 ) 


Similarly for anti-time ordered F W ( t,x ) and 

x), obtained by 6{t) fa 0(—t) in Eo. dBSI) and 
Eo. dBGI) . are 


F^—(t,x) 


F (2) — (t,x) 




x 2 + (a — iv\t\) 2 
a 2 

a — iv\t\ + isgn[t]x 
a — iv\t\ — isgn[t]x 


(B7) 

(B 8 ) 


We absorb the effect of Klein fac¬ 
tor ~i(T c ri R /L{Ti)r] R /L(T 2 )) by introduc¬ 
ing j?( 2 )++/— x) = F( 2 )++/ (t,x) + 

p( 2 )++/ (—t,— x) ± sgn[t]i and f( 2 )h / *~(t,x) = 

2 )+-/ h (t,a:) + F( 2 )h / f(— t,— x) ± i. The general 
form of GW (lu, x) at zero temperature is 


(<*>, x) = 


a isgn[j]k F x 

2 TTCL 


(B9) 


dte^-G* e f ((- ff+ '^)' F<1> ( t ,x)+sgn(j)F ( - 2 ' 1 (t,x)) 


The action 


—Sn = 


dr f dxj [iV9(i, t)8 t <&(x, t) 


-|(K(V0) 2 + -|(V$) 2 )] + - ea)7.} 


At zero temperature 
1 


-1 


^ ($(n)$(r 2 )) = — In 


x 2 + (a + ivt) 


21 


= FW-+(t,x) 
($(ri)0(r 2 )) = ^ In 
= .F ( 2 ) “ + (f,x) 


a + ivt — ix 


a + ivt + ix 


(Bl) 


(B2) 


Here t = t\— t 2 and x = aq—x 2 . For f 2 on the bottom and 
ti on the top branch of Keldysh contour we substitute 
x -A —x and t\ fa t 2 to get 


F (1)+ (t, x) = —7 In 
FM+~(t,x) = -±- In 


x 2 + (a — ivt) 2 


a — ivt + ix 
ao — ivt — ix 


(B3) 

(B4) 


It is straightforward to show that the ’ + (t,x) ob¬ 
tained in Eo. dBDI) (before performing the Fourier trans¬ 
form to frequency space) is the same as Eq.(Al) in 
Ref. USS. 


G ^(t,x) = 


a isgn[j]k F x 

2na 


sgn\j]x — v(t ± * 0 +) 


sgn[j]x + v(t ± * 0 +) 
K+I 


(B10) 


G (t,x) = 


0 isgn[j]k F x 

2na 

a 


sgn[j]x + v(t =F isgn[t] 0 +) 

-i re+i 


sgn[j]x — v{t isgn[t\ 0 +) 


(BH) 


To compute Eci. dBQI) let us first define I-+(u>,x) and 
I~(oj,x) as 


1 1 + 
+|l 

IT 

II 

/•OO 

/ dte iut 


/ —OO 

I-- (w,x) = 

/•OO 

/ dte iut 


J —OO 


x 2 — v 2 (t ± i 0+) 2 
a 2 

x 2 — u 2 (t =p isgn[t\ 0+) 2 


(B12) 
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From Ea. (IB12l) it is easy to check that Ea. (IB9l) is ex¬ 
pressed as 




G (W, *) = 


_ e isgn[j]k F x 


4 7t(k — i) [ a »»"b> + * 


.UJ ~ fl 


e isgn[j]k F x 

4tt(k - §) 


S0n[j]a 


V 

UJ — n 

V 


I-+(u}~ n,x) 


I-- (w - /z,x) 


Since w is real we have (I + (—uj,x))* = I + (uj,x) and 
(/++( — w , x ))* = I (uj,x). We only need to evaluate 
/ ++ and /"* . For I -1 (u>,x) the nonzero contribution 

comes from the lower half circular contour. 



(T 


2 k— i r°° 


FIG. 11: Left: Contour chosen to evaluate I + (ui,x). Right: 
Contour chosen to evaluate I ++ (lo,x) 


dt- 


. M 


— e 


((!) 2 -(* + 2 °+) ; 

e iuy dy 


. — lag I 


—e 


'Cl [~y(y + 2\x\/v)] K l 
e iujy dy 


The full expression for I ++ (uj,x) is 

2aoy/ir f /2i|x|n N 1 K 


T 


= l 


c 2 [y{-y + z\x\/v)\ K 2 

1 —AC 

K K -i(\x\w/iv) 


2(u/tt f /2*|x|u 


zT© — i) \ \ wo 


I ++ (lo,x) = i 
x6(u>) + 


wr 


2 > 
1 —AC 


vY(k — i) I \ WO 2 

—2i\x\v 


K K -i(\x\uj/iv) 


K K -\{—\x\w/iv)6{—uj) >(B16) 


f—2i\x\v \ 1-K 'l 

-!—!_ ) K k _i(— \x\uj/iv) >6(—uj) (B13) We combine the above results and use the derivative 

V w « 2 / J relation 


Here K n (z) is the modified Bessel function of the sec¬ 
ond kind and T(x) is the Gamma function. For evalua¬ 
tion of 7 ++ (w,x), notice that I ++ (uj,x ) is an even func¬ 
tion of uj following its definition: 

I ++ (lo , x) 


( x \—^ / x \—^ 

-J K n (ax) = —a K n+1 (ax) 


(?) 


2 k-i r°° 


dt- 


((f) 2 -(t-* S5 n[t]0+) 2 )' 


(B14) f° r m °dified Bessel functions. Replacing w by w - /i 
to account for nonzero chemical potential and after some 
_ algebras we get 


Thus we only need to evaluate uj > 0 in I ++ (oj,x). For 
this uj > 0 region we have 

I ++ (uj , x)6(ui) 

2 k— i ( r°° / — o \ 


G 4,j ( UJ ’ X ) = 


0 isgn[j]kpx 


i(uj — / Li ) I / 2*| 


2^/HY (k + 7) u 2 1 \(uj - y)a 2 


1 —AC 


(?)' 

r 0 


dte 1 ' 


x \(uj - /z) 


(?) 

K +l. 

|6»(w) 


+ J dte iut - (t + *0+) 

+ J™ dte~ iuit ((f)"-(*“»' 0+ ) 2 ) +2 } 0 H 


x AT, 


— K+i 


Kk-1 

—2i\x\v 
(w — /z)a 2 

*1© - m) 

—iv 


IV 

1 —AC r 


sgn[jx]A" K 


1*1© - /z) 




AT k _i 


>0(/z — w) 


-m) 


- sgn[jx] 


(B17) 


G J"©>*) = 


— e isgn[j]kpx a ( w _ jU ) J / 2*|x|z; 


/ 

.7 cx 


— e 


(?) { 

ie~ wv dy 


.M. 


e iu>v dy 


2^/ttT(k+\) v 2 ) \(uj - fj,)a 2 


(y) 2 + v 2 


'c 3 [~y{y + 2\x\/v)} K 2 
/■°° 


JO 


(y) 2 + v 2 


}d(w) 


Kk-1 


1*1 (W - d) 


= l 


2a-©7r /2*|©*A 77 K _ 1 (|x|w/w)6»(a;) (B15) xi ^- 

nr(«; — 4) V wa / 


—2*|x|u 
(uj — y)a 2 

x\(uj - /z) 


1 —AC 


sgn[3'x]AT re 


1*1 (w - m) 




K k -i 


1*1 (w - /z) 


- S5n[jx] 


>d(cu — /z) 


(BIB) 
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e i S gn[i]fc Fa; a ( w - n) f { 2i\x\v 


2y/nT (k + i) w 2 


(w — /i)a 2 


1 — /■£ 


('tfcil') +sgn \j x]K , fMfcpii> 

\ iv J \ iv 

x9(uj — n) + (same expressions with |x| —> —\x\ 


but lattice constant a oc 1/A). Thus these diver¬ 
gence terms can be safely neglected or suppressed by 
the regularization via a further differentiation^. Us¬ 
ing further the functional relation of Gamma function 
T(2z) = (27r) -1 / 2 2 2z- 3r(z)r(2: + i) we can get Eq. (fl2l) 
from Ea. (lB17|) to Eo. (lB20|) which is consistent with the 
direct derivation done in Ref. 36. 


and sgn[jx\ —> sgn\—jx\)9(g — to) 


(B19) 


Appendix C: Derivation for Green functions 


G to ( w > *) = 


_ e isgn\j]k F x a ( w _ n) j / 2i\x\v 


2^Fr (« + i) 


( u> — g)a 5 


(tfcA) + .gnMK, 

\ IV J \ IV 

x9(g — ui) + (same expressions with |x| -A —|x| 

and sgn[jx\ -A sgn[—jx])9(u} — g) > 


Here we use the Langreth rule^S on the Eo. (fl5l) and 
Eo. dTUl) . The retarded and advanced Green’s functions 
are decoupled from lesser and greater ones and are solved 
directly from these two coupled equations. The explicit 
expressions of all retarded Green’s functions expressed 
via unperturbed Majorana Green’s functions and lead 
Luttinger Green’s functions are: 


(B20) 


K <*{z) - 2 


r(a)(-)“ + T(-a)(-r 

z z 


(1 + 0(s 2 )XB21) 


for z —> 0. For a > 0 the hrst term in Eo. (IB21|) is di¬ 
vergent, reflecting the artifact of the lack of small dis¬ 
tance cutoff in taking the continuous limit (i.e. the 
smallest distance we can take should not be x = 0 


/~i(0)R _ I j. |2 /'lit (/^i{0)R^(0)R _ ^(0)R y~i(Q)R \ 

R _ '-'la \ L P\ ljr 7 P LT 7c73'- T 7/37al 


The x — > 0 limit is obtained by noting that the small 
argument expansion of the modified Bessel function of 
the second kind K a (z) takes the following form^: 


= 


s^iR 
^ 7 a 7/3 


/^iR _ 

W/3 _ 


s^iR 
^7/3 7a 


fR 
J 7 


(Cl) 


/^(0 )R I j. j.* /~iR \ 

^7a 7/3 I V U '7a ^7/3 Lz 7a7/3 Lr 7/37a J 


_ U |2 riR \ 

^7/3 \ L CL | 7a ^7/3 °7a7/3 ^ 7/3 7 a / 


J 7 


(C2) 


^y(0)i? _i_ J. \ 

LT 7^7a “I" i a L /3'^ r 'lp8'lp a V LT 7a ^7/3 Lr 7a7/3 Lr 7/3 7a J 


f? 


Here the numerator //* is defined as 


J 


rR — -i r' i(0)R\j. |2 s~iR j. j.* /^i(0)R s^iR Ij. |2 s~i(0)R/^iR ± /^i(0)R /^iR 

/ 7 = 1 — (_t 7q \l a \ “^a Lr 7 a7/3 Lr ^/3^a "" 1^1 _ ^^S^c^aV'/S 


v //^(0)i? s~i(0)R s^i(0)R /^(0)R \(/^iR s~iR s-^R s~iR \ 
X \ 7a ^7/3 ^ 1 ol 10 ^ 1(^10 ^^a ^0 ^ ^ 0 ^a / 


|ta| 2 |*/»| 2 


(C3) 


Once we obtain the full advanced and retarded Green’s 
functions we then substitute these expressions into the 
equations for lesser and greater ones, which are coupled 


with advanced and retarded Green’s functions. The full 
expressions for lesser Green’s functions are 
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= 


Nll a 1 = 

X 

+ 


Nur 


Nu a 2 = 


De i = 


De 2 = 


Nu a i _ 

De\De 2 De 2 

(1 - t a t}G^ p G^l - \tp\ 2 G$ p G <® R )((1 - \t a \*G$ a G<» R - tltpG^G^l) 

{G®< +G*[\t a \*(GW<G$ a +GM*G<J + tMGM<G^ a +G^ 

(|i a | 2 G£,GW£ + t* a tf,G$, ^G^ r ){G^< + G^[\t a \ 2 (G^<G^ + G^ R G<J 
ft(G' 0 i<G^+G^G<^)] + G; 7 J l^| 2 (G^<G^ +GW*G^) 

t at*p{G^ < G^ a ^ f) + G^ fl G| a ^)]}^ 
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